Casimir interactions in graphene systems 
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The non-retarded Casimir interaction (van der Waals interaction) between two free standing 
graphene sheets as well as between a graphene sheet and a substrate is determined. An exact 
analytical expression is given for the dielectric function of graphene along the imaginary frequency 
axis within the random phase approximation for arbitrary frequency, wave vector, and doping. 
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The first reference to the material graphene in the liter- 
ature was made by Boehm et. al[U in 1962. With modern 
technology it is possible to produce large area graphene 
sheets and graphene has become one of the most ad- 
vanced two-dimensional (2D) materials of today. It was 
awarded the nobel prize in physics in 2010. Due to its 
superior transport properties it has a high potential for 
technological applications!!^. A free standing graphene 
sheet has a very interesting band structure. The valence 
and conduction bands form two sets of cones. In each 
set the two cones are aligned above each other with their 
tips coinciding at the fermi level. Thus, the fermi sur- 
face is just two points in the Brillouin zone; the value 
of the band gap is zero. The energy dispersion in the 
conduction and valence bands is linear which means that 
the carriers behave as relativistic particles with zero rest 
mass. When a graphene layer is formed on a substrate 
the fermi level moves up or down in energy, the sheet is 
doped. A graphene layer interacts with other graphene 
layers or with a substrate with CasimirQ forces. These 
forces are very important in graphene systems. Many- 
body effects also modifies the dispersion of the bands|8- 
Il3j |. The present work is devoted to the forces. 

We begin by calculating the interaction energy between 
two graphene layers. We limit the calculation to small 
separations where retardation effects are negligible. Thus 
we calculate the non-retarded Casimir interaction or in 
other words the van der Waals interaction. The interac- 
tion energy is then nothing but the inter-layer correlation 
energy [14]. At zero temperature it is given by 





Figure 1. Feynman diagrams for the correlation energy in the 
two graphene sheet system. The ellipses represent polariza- 
tion bubbles and the dashed lines the interactions indicated 
in the figure. The numbers 1 and 2 refer to which sheet the 
electron belongs to. See [TJ for details. 



feynman diagrams representing this energy are given in 
Fig. [1] The force is obtained as minus the derivative of 
the energy with respect to separation, d, i.e., 
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where cvo (<?, uj) is the polarizability of one graphene layer. 
The prime indicates that the function is calculated along 
the imaginary axis of the complex frequency plane. In 
terms of the polarizability the dielectric function is given 
by £ (q,ui) = l + a (q,u) = 1 - v 2D (q) x (q, w)/k, where 
v 2D (q) = 2ire 2 j q is the two-dimensional fourier trans- 
form of the coulomb potential, k is the dielectric constant 
of the surrounding medium, and \ (<Zi w ) the density- 
density correlation function or polarization bubble. The 
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The retarded version of the density-density correla- 
tion function, x(g,w), for graphene has independently 
been obtained by Hwang and Das Sarma[15j and by 
>, Wunsch[16] et al. Here we present our general expres- 
sion in the complex frequency plane, away from the real 
axis. In particular we give the result on the imaginary 
axis, which is where we perform the present calculations. 
We, just as in include the contribution from 

the conduction and valence bands, only, and assume that 
the linear dispersion of the bands extends forever. In the 
real system the dispersion starts to deviate from linear at 
some point and transitions from the occupied core states 
to empty bands higher up in energy will contribute to 
some extent to the screening. We neglect this here. Let 
us first begin with an undoped graphene sheet. In a gen- 
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eral point in the complex frequency plane, away from the 
real axis the density-density correlation function is 



x(q>z) = - 



WH y / v 2 q 2 
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where v is the carrier velocity which is a constant in 
grapene (E — ±hvk), and g represents the degeneracy 
parameter with the value of 4 (a factor of 2 for spin and 
a factor of 2 for the cone degeneracy.) 

With this particular screening it turns out that 
a' (q/<i, uj/d) — a (q, lu) and the separation dependence 
of the interaction becomes very simple, after a change of 
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10 5 to/s; 



with the value for v chosenfl6| to be 8.73723 
d is the distance in A. 

When the graphene sheet is doped the expression 
for the density-density correlation function is much 
more complicated. In the two next equations we use 
dimension-less variables. The variable x is the momen- 
tum divided by 2k f] the variable y is the frequency di- 
vided by 2Ep/h; the variable z is a general complex val- 
ued frequency divided by 2Ep/h. 

In a general point in the complex frequency plane, 
away from the real axis the density-density correlation 
function is 
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On the imaginary axis it is 

x' (q, £*>) = x: (q> iw ) = -Do { 1 + 



(5) 



4yV 



7T — atan 



j{ [^(^-l) + (^+l) 2 ] 2 +(2^ 2 ) 2 } sin{i 



2 (y 2 -i) + (v 2 +W 



v / {x 2 +y 2_ 1) 2 + {2y) 2^ (y 2 + 1) 



\f^2x 2 (y 2 -l)-2(y i -6y 2 + l)+2(y 2 + l)y/x i +2x 2 (y 2 -l) + (y 2 + l^ 



where the density of states at the fermi level, Do, is 
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Here n is the doping concentration. The same result 
holds for excess of electrons and excess of holes. 

The numerical results for the size of the interaction 
energy between two undoped graphene layers in vacuum 
is shown as the straight line in Fig. [5J The bent curves 
are valid for doping concentrations 1 x 10 10 , 1 x 10 11 , 1 x 
10 12 ,andl x 10 13 cm~ 2 , respectively, counted from be- 
low. The interaction energy is negative, leading to an 
attractive force. The interaction between two 2D metal- 
lic sheets was derived in [14]. To illustrate the difference 
between the two systems we show in Fig. [3] the corre- 
sponding results for two 2D metallic sheets in vacuum 
for the same set of carrier concentrations; we have used 
the effective mass unity for the carriers. The van der 
Waals interaction between two atoms varies with sepa- 
ration, d, as d~ e . Between two half spaces it varies as 
d~ 2 . We found in [TiJ that the system of two 2D metal- 
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lie sheets has a separation dependence with a fractional 
power, v.i.z., <i~ 5 / 2 (see the straight part of the curves in 
Fig. [3]) Here, we have found a system with yet another 
separation dependence. For two undoped graphene layers 
the van der Waals interaction varies as d~ 3 . This is the 
same dependence as for the Casimir interaction between 
two half spaces. The origin of the half-integer behavior 
for the 2D metal sheets is the square root behavior of the 
dispersion curve for 2D plasmons. In undoped graphene 
the plasmon dispersion curve is linear. In the doped case, 
however, the plasmon dispersion curve attains the square 
root behavior and for larger separations the interaction 
varies as <i~ 5 / 2 (see the rightmost part of Fig. For 
small enough separations the contribution to the interac- 
tion from the free carriers saturates. This is obvious in 
both Figs. HJandEl 

Before we treat next geometry which is a graphene 
sheet above a substrate it is illustrating to rederive Eq. 
(fTJ) in a different way. In [3| we derived it in two al- 
ternative ways; here we do it in yet another way. Let 
us assume that we have an induced carrier distribution. 
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Figure 2. The attractive interaction energy between two 
graphene sheets. The straight line is for undoped sheets, while 
the bent curves are for doping densities 1 x 10 10 , 1 x 10 11 , 1 x 
10 12 ,andl x 10 13 cm , respectively, counted from below. 
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Figure 3. The attractive interaction energy between two 
2D metallic layers. The curves are for doping densities 
1 x 10 1 

counted from below. 



1 x 10 n ,l x 10 12 ,andl x 10 13 cm" 2 , respectively, 



Pi (q, lu), in sheet 1. 
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This gives rise to the potential 



v (^q, lu) — v""' (q) pi (q, a;) in sheet 1. It results in the 
potential exp (— qd) v 2D (q) p\ (q, lu) in sheet 2, which is 
parallel with sheet 1, the distance d away. This potential 
is screened by the carriers in sheet 2. The resulting po- 
tential is exp (—qd) v 2D (q) p\ (q, lu)/[1 + cto (q, lu)]. This 
gives rise to an induced carrier distribution in sheet 2, 

P2 (q, U ) = X (q, lu) e-« d v 2D (q) g ^ . (8) 

1 + a (q,w) 

In complete analogy, this carrier distribution in sheet 2 



gives rise to a carrier distribution in sheet 1, 

Pl (q, U ) = X (q, «) e"* V* (q) P ^f . (9) 

1 + a (q,w) 

To find the condition for self-sustained fields, normal 
modes, we let this induced carrier density in sheet num- 
ber 1 be the carrier density we started from. This leads 
to 
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The left hand side of this equation is exactly the argu- 
ment of the logarithm in Eq. ([I]). This equation was 
derived using many-body theory. In [t7| the interaction 
energy is derived from the electromagnetic normal modes 
of the system. One ends up with an identical expression 
to the one in Eq. ([1} where now the argument of the log- 
arithm is the function in the condition for normal modes. 
If we have a 2D layer (like a graphene sheet) the distance 
d above a perfect metal substrate the procedure is very 
similar. We start with an induced mirror carrier density, 
pi (q, lu), in the substrate. The induced carrier density in 
the graphene sheet is given by the expression in Eq. f8j) 
except that now the distance beween the mirror charge 
and the graphene layer is 2d instead of d. Eq. @ is then 
replaced by p\ (q, lu) = —p 2 (q, lu), according to the re- 
sult for a mirror charge at a perfect metal interface. The 
condition for modes becomes 
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Unfortunately, for a graphene layer above a perfect metal 
substrate the integral does not converge, without a cutoff; 
the polarizability dies off too slowly with frequency; for 
a 2D metallic sheet, on the other hand, above a perfect 
metal substrate there is no problem. 

Next we focus on a graphene layer above a real sub- 
strate. The derivation is the same as for a perfect metal 
substrate until the last step. Now, using the theory of 
image charges we realize that the relation between the 
induced carrier densities is 



1 



Pi (q,^j 



-p2 (q,w) 



and the condition for normal modes becomes 
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We have neglected spatial dispersion in the substrate; in- 
clusion of spatial dispersion would lead to a higher order 
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of complexity [18] and would have negligible effects on the 



present results. From the condition in Eq. (|14[) follows 
that the interaction energy is given by 
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Figure 4. The attractive interaction energy between a 
graphene sheet and a gold substrate. The lowest curve is for 
an undoped sheet, while the other curves are for doping densi- 
ties 1 x 10 1() , 1 x 10 u , 1 x 10 12 , andl x 10 13 cm" 2 , respectively, 
counted from below. 

The result for a graphene sheet above a gold substrate 
is shown in Fig. 2] The dielectric function of gold along 
the imaginary frequency axis was obtained from experi- 
mental data extrapolated in a way described in (l9| and 
using a modified Kramers Kronig dispersion relation (see 
Eq. (6.75) in 0.) 

In summary, we have derived and calculated the 
non-retarded Casimir interaction between two graphene 
sheets and between a graphene sheet and a substrate. We 
found the interaction energy between two virgin graphene 
sheets varies with separation, d, as d~ 3 (the force as d~ A .) 
This is an unexpected result. For a pair of films summa- 
tion over pair interactions leads to a d~ 4 dependence. 
The interaction between doped graphene sheets has a 
more complex separation dependence. We have derived 
the dielectric function of graphene along the imaginary 
frequency axis within the random phase approximation 
for arbitrary frequency, wave vector, and doping. These 
results are needed for the present calculations and for 



future calculations of many-body effects in graphene. 
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